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In this paper I investigate the usability of the characteristic functions for the description of the 
dynamics of open quantum systems focussing on non-Lindblad-type master equations. I consider, as 
an example, a non-Markovian generalized master equation containing a memory kernel which may 
lead to nonphysical time evolutions characterized by negative vàlues of the density matrix diagonal 
elements [S.M. Barnett and S. Stenholm, Phys. Rev. A 64, 033808 (2001)]. The main result of the 
paper is to demonstrate that there exist situations in which the symmetrically ordered characteristic 
function is perfectly well defined while the corresponding density matrix loses positivity. Therefore 
nonphysical situations may not show up in the characteristic function. As a consequence, the 
characteristic function cannot be considered an alternative complete description of the non-Lindblad 
dynamics. 
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The theory of open quantum systems describes the in- 
teraction of a quantum system with its environment 0. 
Although many physical systems may be considered un- 
der certain conditions quasi-closed for certain intervals 
of time, no quantum system may be seen as completely 
isolated from its surroundings. The unavoidable interac- 
tion between the system and its environment leads to the 
phenomena of decoherence and dissipation . 

The description of the dynamics of open systems has 
recently attracted much of attention for mainly two rea- 
sons. On the one hand environment induced decoher- 
ence due to the establishment of entanglement between 
the system and the environment is one of the key issues 
of the quantum measurement theory 0. On the other 
hand, the system-environment interaction seems to be 
the major limiting factor in the realization of quantum 
devices necessary for the new quantum technologies, e.g. 
quantum computation |3J. 

The study of the dynamics of an open system is, in 
general, a very difhcult task. Usually, even if one is in- 
terested in the dynamics of the system only, the influence 
of the typically infinite environmcntal degrees of freedom, 
makes it impossible to solve exactly the equations of mo- 
tion for the observables of interest. For this reason, the 
Standard description of open systems rely on a number of 
approximations which allow to derive a master equation 
for the reduced density matrix of the system. The two 
most common approximations are the weak coupling ap- 
proximation, vàlid when the interaction between system 
and environment is sufRciently weak, and the Markov 
approximation, relying on the assumption that the char- 
acteristic times of the system are much larger than those 
of the environment Q| . 

Generally, when these approximations are satisficd, the 
master equation for the reduced density matrix may be 



written in the so-called Lindblad form, which is the only 
possible form of first-order linear differential equation, 
for a completely positive dynamical semigroup having 
bounded generator 5, 6]. The Lindblad master equation, 
however, is vàlid as long as the weak coupling and Markov 
approximation hold. While these assumptions are of- 
ten well justified in quantum òptics, in many solid-state 
systems, i.e. photonic band-gap materials and quantum 
dots, the Markov approximation does not hold [?|. Sim- 
ilarly, the reservoir interacting with a single mode cav- 
ity in atom làsers is strongly non-Markovian ||. Non- 
Markovian generalized master equations usually ar not 
of Lindblad type. 

It is worth noticing that there exist also Markovian 
systems described by master equations which cannot be 
cast in the Lindblad form 0, llfll ] . An important problem 
in the description of open quantum systems whose mas- 
ter equations are not in the Lindblad form is that their 
dynamical map needs not be completely positive 
and this may lead to physical inconsisteney. It may even 
happen that the positivity condition of the density ma- 
trix during the time evolution, a condition less restrictive 
then complete positivity but necessary to guarantee the 
probabilistic interpretation of the density matrix, breaks 
down. 

For the sake of completeness, let me underline that 
complete positivity is a necessary requirement for a con- 
sistent physical description of open quantum systems 
whenever factorized initial conditions for the system and 
the reservoir are assumed, i.e. pr(0) = /3(0)<8>/3e(0), with 
Pt(0), /5(0), j5e(0), initial density matrices of the total sys- 
tem, of the reduced system of interest and of the envi- 
ronment, respectively. Most of the derivations of master 
equations found in the literature rely on this assump- 
tion. However, when correlations are present at the ini- 
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tial time, acceptable quantum dynamics which are not 
completely positive may exist |l3j . In the following I will 
focus on the case of factorized initial conditions for which 
the dynamical map must be completely positive. 

When working with non-Markovian generalized mas- 
ter equations, or with master equations which are not in 
the Lindblad form, it is of crucial importance to establish 
conditions under which the density matrix preserves pos- 
itivity and complete positivity during the time evolution. 
In most of the cases these conditions are given in terms 
of the density matrix elements at time t, and therefore 
rcquirc the knowledge of the analytic solution of the non- 
Markovian master equation. In Ref. [ï^. conditions for 
complete positivity in terms of the memory kernel are 
presented for a class of generalized master equations. To 
the best of the author's knowledge, however, a criterion 
analogous to the Lindblad one for general master equa- 
tions with memory has not yet been formulated. 

Very useful tools for the description of the dynam- 
ics of paradigmàtic open quantum systems such as the 
damped harmònic oscillator (quantized mode of the elec- 
tromagnètic field, motion of a trapped ion) or the quan- 
tum Brownian particle, are the characteristic functions 
and the quasi-probability distribution functions. Both of 
them contain all the information necessary to reconstruct 
the density matrix, and therefore they have been consid- 
ered up to now " alternative complete descriptions of the 
dynamics "0,0- However, the characteristic function 
can be considered an alternative complete description of 
the dynamics if and only if it is equivalent to the den- 
sity matrix. I will prové in the following that, contrar- 
ily to what has been believed until now, the character- 
istic function and the density matrix descriptions of the 
dynamics cannot be considered equivalent. To the best 
of the author's knowledge this is the first time that the 
equivalence between the characteristic function and the 
density matrix is questioned, and an example pointing 
out the non-equivalence between these two approaches is 
presented. 

Let me begin by recalling the definition of the p- 
ordered characteristic functions 



In what follows I will focus on the p = characteris- 
tic function x(£iP — 0) = x(£)' known as simmetrically 
ordered characteristic function (SCF) or quantum char- 
acteristic function. Having in mind Eq. Q it is straight- 
forward to prové that the result obtained in this paper 
for the simmetrically ordered characteristic function also 
applies to the other two characteristic functions. 

The SCF is always defined and it is, in general, a 
complex- valued function satisfying the following proper- 
ties: 



X(É = 0) = 1; 



\x(0\ < i- 



(2) 



The first of the two properties is a consequence of the 
fact that Tr[/5] = 1, while the second stems from the fact 
that x(£) is the expectation value of the displacement 
operator D(Ç) which is unitary, and therefore the magni- 
tude of its eigenvalues is unity. One of the advantages of 
using the simmetrically ordered characteristic function is 
that the analytic expression for the mean vàlues of many 
observables of interest may be calculated easily by means 
of the relation 



(3) 



In the literature the characteristic functions have been 
extensively used to study the dynamics of both Marko- 
vian [iMllaLÜÏ and non-Markovian [0$ Lüï \M uÈ uÈ 
open systems. It has also been shown that characteristic 
functions may be used to establish observable conditions 
of nonclassicality for the states of the quantized electro- 
magnètic field [23j . In what follows I will show that there 
exist situations in which the characteristic function de- 
scription of an open quantum system may lead to prob- 
lems. 

Let me consider the non-Markovian dynamics of a 
harmònic oscillator interacting with a zero temperature 
reservoir. I consider one of the most popular phenomeno- 
lo gica l model for this systems, involving a memory kernel 
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X(f,p) = Tr [pD(0\ exp (p\Ç\ 2 /2) 

= Tr[pexp(£àt-ra)]exp(p|£|72), (l) 

where D(Ç) is the Glauber displacement operator and 
à (éft) is the annihilation (creation) operator of the 
quantum harmònic oscillator. In the previous equa- 
tion, the parameter p assumes the vàlues p = 1,0,-1 
in correspondence to normal, symmetric and antinor- 
mal ordering of the creation and annihilation operators. 
The two-dimensional Fourier transform of xÍÇjP) gives 
the Glauber-Sudarshan P-representation for p = 1, the 
Wigner function for p — 0, and the Husimi Q-function 
forp= -1 \Ü&. 



dp(t) 
dt 



K(t - t')Cp{t')dt', 



(4) 



where K(t — t') is the memory kernel and the Liouvillian 
operator C is given by 



Cp = làpa} 



à'àp — pa'à. 



(5) 



By applying the rules 



- l d t\ "t- ( d ?\ 



pa 



de 



dC 



(6) 
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one may derive, from the master equation Q, the corre- 
sponding integro-difFerential equation for x(£) 



dt 



I K(t - t 1 ) x 



Following a method developed in |2j| for non-Markovian 
Fokker-Plank equations, a formal solution of this integro- 
difFerential equation may be obtained in form of an in- 
tegral decomposition involving the solution of the cor- 
responding Markovian problem, which is known in the 
literature (see, e.g., [ïjjp. 

Let me focus on the case of a memory kernel of expo- 
nential type 



K{t-t')=g 2 e-^ t - t '\ 



(8) 



with g coupling strength and 7 decay constant of the 
system-reservoir corrclations. I consider, as initial state, 
a Fock state \n), with \n) being the eigenstates of the 
quantum number operator h = à'à. The corresponding 
SCF reads as follows 



Xn(0 = L n m 2 )e 



-|í| 2 /2 



(9) 



with L„(|£| 2 ) the Laguerre polynomial of order n. For 
n = 1, e.g., the initial simmctrically ordered characteris- 
tic function is x(£,0) = (1 — |C| 2 ) e ^' an d it is easy 
to verify by direct substitution that 



- 7 i/2 



2Ü 



sin Í7í 



,-|í| 2 /2 



(10) 

is a solution of Eq. 0) , with SI = \j2g 2 — (j/2) 2 . Figure 
n shows the absolute value of the simmetrically ordered 
characteristic function, as given by Eq. IjlOfl. as a func- 
tion of |£| 2 at different time instants and for g/7 = 1. 
Looking at Eq. í|lf)|> it is easy to verify that, for the pa- 
rameters considered in the example, the SCF is always 
dcfincd and it satisfies at all times and for each vale of 
Ç the conditions given by Eqs. J2J, as one can also see 
clearly from the figure. 

Once the time evolution of the simmetrically ordered 
characteristic function is known one can always recon- 
struct the density matrix at all times, since p(t) can be 
obtained from x(£>^) by using the relation 



1 

2^ 



(11) 



From this equation, keeping in mind that the diagonal 
elements of the Glauber displaccment operator are given 
by 




(n\D(0\n)^L n (\e)e-^ 2/2 , 



(12) 



FIG. 1: Behavior of |x(|Ç| 2 ) | , as given by Eq. (TÜJ, at different 
times instants in the interval < r < 0. The initial state is 
\n — 1} and 3/7 = 1. 



and that ii(|Ç| 2 ) = 1 - |£| 2 , one gets 

p n (t) = (n = l|p(í)(n = 1) 

= e~ 7í/2 (cos fit + ^ sin Qt^j . (13) 

The quantity p\\(t) describes the decay of the popula- 
tion of the initial state \n = 1) due to the interaction 
with a zero-temperature reservoir. For t — > 00, p\i — > 
and /Soo — * 1, indicating that, due to dissipation, the 
state of the oscillator passes from the initial excited Fock 
state \n = 1) to the ground state. A close look to Eq. 
(ÇHfll . however, shows that the positivity of the density 
matrix is clearly violated since pu(t) < for some inter- 
vals of time. This situation has been discussed in detail 
by Barnett and Stenholm in [ï^, where the risks of an 
apparently physically well grounded memory kernel, as 
the exponential one given by Eq. JHJ), were carefully an- 
alyzed. The loss of positivity shows up for strong or 
intermediate coupling regime ng 2 /~f > 1/8 [l2] | alerting 
us about the fact that the dynamics of the system, for 
these vàlues of the parameters, is unphysical. 

It is worth stressing, however, that if one describes the 
time evolution by using the simmetrically ordered char- 
acteristic function, one does not realize that the equa- 
tions of motion lose physical sense because this function, 
contrarily to the density matrix which violates one of its 
defining conditions (positivity), continues to verify at all 
times the conditions given by Eqs. @. Hence for non- 
Lindblad cases, when both the complete positivity and 
the positivity conditions may be violated, unphysical sit- 
uations, such as the negativity of the density matrix, may 
not show up in the dynamics of the SCF. 

A careful analysis of Eq. I|1U|) shows that only for vàl- 
ues of the ratio 3/7 such that, for certain intervals of 
time, pn(í) < —1/2 then the second of the conditions 
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given in Eqs. © is violated for |£| 2 <C 1. In general, 
however, there is no correspondence between the loss of 
the positivity condition, and therefore of complete posi- 
tivity, and the violation of one of the conditions defming 
the simmetrically ordered characteristic function. Stated 
another way the problem is the following. The fact that 
£>(£) is an unitary operator, together with the properties 
that p is a positive (thcn Hcrmitian) trace-class operator 
with trace 1, imply that x(£)l < 1> l - c - the second of Eq. 

This is, however, only a necessary condition, indeed 
I have shown in the paper that there exist situations in 
correspondence of which the density matrix is not posi- 
tive but still |x(£)| < 1. The crucial question is therefore, 
which is the additional condition to be imposed on the 
SCF to ensure that the operator p, defined through Eq. 
(|llf). is a positive trace-class operator with trace 1? The 
answer to this question is very important since it would 
allow to use safely the SCF in non-Lindblad cases. Un- 
til when this condition is not found, one cannot claim 
that the density matrix and the SCF are equivalent de- 
scriptions of the dynamics, as the example given in this 
paper clearly indicates. The derivation of the condition 
to be imposed on the simmetrically ordered characteris- 
tic function to make it a useful tool in the description 
of non-Lindblad type situations will be the object of fur- 
ther study. However, presently, it seems to the author 
that this question does not have a simple answer. 

Let me conclude considering the behaviour of the 
Wigner function. The Wigner function is the two- 
dimensional Fourier transform of the simmetrically or- 
dered characteristic function 



W(a) = \ 



aÇ* — a* Ç 



(14) 



The Wigner function is a real valued function satisfying 
the condition |W(a)| < 2/tt. Inserting Eq. l(ÏÜ|l into Eq. 
(tHjl we get 



W(a,t) = -e 

7T 



[1 + 2(2|a| 2 - l)pn(í)] , (15) 



with pu(t) given by Eq. I|13fl . From the previous equa- 
tion one may verify that whenever < pn{t) < 1, then 
W^oOI < 2/7r, but the former inequality is immediately 
violated in correspondence to a violation of the positivity 
condition, in our example when pn{t) becomes negative. 

It is worth stressing the difference between the Wigner 
function, which is the Fourier transform of the SCF, and 
the simmetrically ordered characteristic function itself. 
While the first one violates one of its defming conditions 
when the density matrix loses positivity, the second one 
does not. In this sense it seems that the simmetrically 
ordered characteristic function has less "physical mean- 
ing"than the density matrix or the Wigner function. It 
is worth noticing that, as a consequence of the lack of a 
condition equivalent to the positivity of the density ma- 
trix, each time one deals with non-Lindblad dynamics 



one can use the SCF only if it is possible to derive the 
corresponding density matrix, by means of Eq. I|ll|) , and 
check its positivity. In this paper I considered a rather 
easy example of the dynamics for which both the den- 
sity matrix and the SCF solutions have simple analytic 
expressions. In general, however, it is not obvious that, 
once the solution of the Fokker-Planck equation for the 
SCF is known, one is also able to derive a useful ex- 
pression for the density matrix necessary to check the 
positivity condition. This fact strongly limits the usabil- 
ity of the SCF for the study of non-Lindblad dynamics. 
For this reason new necessary and sumcient conditions 
establishing the cquivalcncc between the SCF and the 
density matrix are highly desirable. Moreover, the an- 
swer to the open question posed in this paper about the 
idcntification of a condition on the SCF correspondent to 
positivity of the density matrix would shed light on the 
physically meaningful ingredient allowing to consider the 
description of an open quantum system via the SCF as a 
complete description of its dynamics. 

This work has been supported by the European 
Union's Transfer of Knowledge project CAMEL (Grant 
No. MTKD-CT-2004-014427) and by the Italian Na- 
tional Foundation Angelo Delia Riccia. The author 
gratefully acknowledges Antonino Messina, Francesco 
Petruccione, Jyrki Piilo, and Stig Stenholm for the useful 
comments and suggestions about the paper. 



[1 
[2 

[3 

[4 

[5] 

[6] 

[<: 

[§ 
[9 
[10 
[11 
[12 
[13 



Electronic address: maniscalco@ukzn.ac.za 

H. -P. Breuer and F. Petruccione, The Theory of Open 
Quantum systems (Oxford University Press, 2002). 

E. Joos, H.-D. Zeh, C. Kiefer, D. Giulini, J. Kupsch, and 

I. -O. Stamatescu, Decoherence and the Appearence of a 
Classical World in Quantum Theory, 2nd ed. (Springer- 
Verlag, Berlin, 2003). 

M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, England, 2000). 

R. Alicki and K. Lendi, Quantum Dynamical Semigroups 
and Applications (Springer-Verlag, Berlin, 1987). 
G. Lindblad, Commun. Math. Phys. 48, 119 (1976). 
V. Gorini, A. Kossakowski, and E.C.G. Sudarshan, J. 
Math. Phys. 17, 821 (1976). 

S. John and T. Quang, Phys. Rev. Lett. 78, 1888 (1997); 
T. Quang, M. Woldeyohannes, S. John, and G. S. Agar- 
wal, Phys. Rev. Lett. 79, 5238 (1997). 
J.J. Hope, GM. Moy, M.J. Collett, and CM. Savage, 
Phys. Rev. A 61, 023603 (2000). 

W.G Unruh and W.H. Zurek, Phys. Rev. D 40, 1071 
(1989). 

W.J. Munro and GW. Gardiner, Phys. Rev. A 53, 2633 
(1986). 

K. Kraus, States, Effects and Operations, Fundamental 
Notions of Quantum Theory (Acadèmic, Berlin, 1983). 
S.M. Barnett and S. Stenholm, Phys. Rev. A 64, 033808 
(2001). 

P. Pechukas, Phys. Rev. Lett. 73, 1060 (1994); R. Alicki, 



■5 



ibid. 75, 3020 (1995); P. Pechukas, ibid. 75, 3021 (1995). 
[14] A.A. Budini, Phys. Rev. A 69, 042107 (2004). 
[15] S.M. Barnett and P.M. Radmore Methods in Theoreti- 

cal Quantum Òptics (Oxford University Press, Oxford, 

1997). 

[16] L. Mandel and E. Wolf, Optical Coherence and Quantum 
Òptics (Cambridge University Press, Cambridge, 1995). 

[17] V. Giovannetti et ai, Phys. Rev. A 70, 032315 (2004). 

[18] A. Isar et al, Int. J. Mod. Phys. E, Vol. 3, No 2, 635 
(1994). 

[19] W.T. Strunz and F. Haake, Phys. Rev. A 67, 022102 
(2003). 



[20] F. Intravaia, S. Maniscalco, and A. Messina, Phys. Rev. 

A 67, 042108 (2003). 
[21] S. Maniscalco, J. Piilo, F. Intravaia, F. Petruccione and 

A. Messina, Phys. Rev. A 69, 052101 (2004). 
[22] S. Maniscalco, J. Piilo, F. Intravaia, F. Petruccione and 

A. Messina, Phys. Rev. A 70, 032113 (2004). 
[23] W. Vogel, Phys. Rev. Lett. 84, 1849 (2000); Th. Richter 

and W. Vogel, ibid. 89, 283601 (2002); E. Shchukin, 

Th. Richter, and W. Vogel, Phys. Rev. A 71, 011802(R) 

(2005). 

[24] I.M. Sokolov, Phys. Rev. E 66, 041101 (2002). 



